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ABSTRACT 

Mars was second to Jupiter in being recognized as the host of a population of Trojan minor 
bodies. Since 1990, five asteroids - 5261 Eureka, (101429) 1998 VF31, (121514) 1999 UJ7, 
2001 DH47 and (31 1999) 2007 NS2 - have been identified as Mars Trojans, one L4 and four 
L5. Dynamical and spectroscopic evidence suggests that some Mars Trojans may be remnants 
of the original planetesimal population that formed in the terrestrial planets region. Here, we 
revisit the long-term dynamical evolution of the previously known Mars Trojans and show that 
201 1 SC191, 201 1 SL25 and 201 1 UNgs are also trailing (L5) Mars Trojans. They appear to be 
as stable as Eureka and may have been Trojans over the age of the Solar system. The fact that 
five Trojans move in similar orbits and one of them is a binary may point to the disruption of a 
larger body early in the history of the Solar system. Such a catastrophic event may also explain 
the apparently strong asymmetry in terms of number of objects (one versus seven) between 
the L4 and L5 regions. Future spectroscopic observations should be able to reject or confirm a 
putative common chemical signature that may lend further support to a collisional scenario. 

Key words: celestial mechanics - minor planets, asteroids: individual: 201 1 SC191 - minor 
planets, asteroids: individual: 2011 SL25 - minor planets, asteroids: individual: 2011 UNga - 
planets and satellites: individual: Mars. 



1 INTRODUCTION 

Trojans are minor bodies that share the semimajor axis of their host 
body but they may have different eccentricity and inclination. In a 
frame of reference rotating with the host, they move in the so-called 
tadpole orbits around the Lagrangian equilateral points L4 and L5 ; 
L4 is located on the orbit of the host object at some 60° ahead or 
east of the host, and L5 is some 60° west. In other words and for 
a Trojan object, the relative mean longitude Ar = A — Ah oscil- 
lates around the values +60° (L4) or -60° (or 300°, L5), where A 
and Ah are the mean longitudes of the Trojan and the host body, 
respectively. The first Trojan object, 588 Achilles, was discovered 
by M. Wolf in 1906 at Jupiter's L4 point (Einarsson 1913). Since 
1990, Mars is second to Jupiter in being recognized as the host of 
a population of Trojan minor bodies (see, for example, Marzari et 
al. 2002 for a review). The list of Mars Trojans currently includes 
5261 Eureka (Bowell et al. 1990; Holt et al. 1990; Mikkola et al. 
1994), (101429) 1998 VF31 (Ticha et al. 1998; Tabachnik & Evans 
1999), (121514) 1999 UJ7 (Rivkinetal. 2003), 2001 DH47 (Scholl, 
Marzari & Tricarico 2005) and (311999) 2007 NS2 (Rodriguez et 
al. 2007; Izidoro, Winter & Tsuchida 2010; Schwarz & Dvorak 
2012). Right after the discovery of Eureka, it was suggested that 
Mars Trojans may have been captured fairly late in Solar system 
history (Bowell et al. 1990) but soon after, the primeval Mars Tro- 
jan scenario started to gain support from both dynamical (long-term 
integrations) and spectroscopic data (Mikkola et al. 1994; Howell 
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et al. 1995). Today, it is widely accepted that some Mars Trojans 
may be primordial bodies, perhaps the only surviving examples of 
the planetesimal population that formed in the teiTestrial planets 
region (e.g., Connors et al. 2005; Scholl et al. 2005). 

Any claim regarding the putative ancient nature of all or part 
of the current Mars Trojan populations must be supported by dy- 
namical stability analyses. The stability of Martian Trojans was 
first studied by Mikkola & Innanen (1994). For a simulated time 
of 100 Myr, they found two regions of stability depending on the 
orbital inclination, i; only orbits with i £ (15, 30)° or i G (32, 44)° 
appeared to be long-term stable. Tabachnik & Evans (1999, 2000) 
found that the inclination ranges in Mikkola & Innanen (1994) en- 
sured stability around the Lagrangian point L4 but that around L5, 
the range of stable inclinations was i G (15, 40)°. In their 100 Myr 
simulations, both Eureka and 101429 were deeply settled into the 
L5 stable region; they also suggested a primordial origin for these 
two objects. Focusing on the role of secular resonances in the dy- 
namical evolution of Mars Trojans, Brasser & Lehto (2002) found 
that several secular resonances with the Earth, Mars and Jupiter 
were able to explain the unstable inclination ranges found by previ- 
ous studies. Scholl et al. (2005) further explored this topic to con- 
clude that secular resonances make Trojan orbits with i < 15° as 
well as those with i ~ 30° unstable. Orbits with i > 35° are also 
unstable due to the Kozai resonance (Kozai 1962). In their 4.5 Gyr 
calculations, they also found that the dynamical half-life of Mars 
Trojans in the most stable regions is of the order of the age of the 
Solar system and that the Yarkovsky effect has negligible role for 
objects larger than about 10 m. All known Mars Trojans have rel- 
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atively small sizes with diameters of order of 1 km (Trilling et al. 
2007). Numerical modelling (Tabachnik & Evans 1999, 2000) pre- 
dicts that the number of Mars Trojans with diameter > 1 km could 
be as high as 30. 

Here, we search for Mars Trojan candidates among known as- 
teroids with relative semimajor axis \a — au3,s\ < 0.001 au, then 
perform A'^-body calculations to confirm or reject their Trojan na- 
ture. In this Letter, we show that the minor bodies 2011 SCigi, 
2011 SL25 and 2011 UNea are stable trailing (L5) Mars Trojans. 
This Letter is organized as follows: in Section 2, we briefly outline 
our numerical model. Section 3 focuses on currently known Mars 
Trojans, reviewing their dynamical status. Section 4 presents our 
results for the objects studied here. These results are discussed and 
our conclusions summarized in Section 5. 



2 NUMERICAL MODEL 

In order to confirm the current Trojan nature of our candidates, 
we perform A'^-body calculations in both directions of time (±2 
Myr) with the Hermite integrator (Makino 1991; Aarseth 2003), 
in a model Solar system which includes the perturbations by the 
eight major planets and treat the Earth and the Moon as two sep- 
arate objects; it also includes the barycentre of the dwarf planet 
Pluto-Charon system and the five largest asteroids, 1 Ceres, 2 Pal- 
las, 4 Vesta, 10 Hygiea and 31 Euphrosyne. Our calculations con- 
sider point, constant mass objects orbiting in a conservative system; 
therefore, relativistic effects are ignored and no modeling of the 
Yarkovsky and Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) 
effects (see, e.g., Bottke et al. 2006) is attempted. The standard ver- 
sion of this direct A'^-body code is publicly available from the loA 
web siteQ additional details can be found in de la Fuente Marcos 
& de la Fuente Marcos (2012). In order to validate our results, we 
have also calculated the evolution of the previously known Mars 
Trojans (see Figs[T]and|2}. Our orbital integrations use initial condi- 
tions (positions and velocities in the barycentre of the Solar system 
referred to the JD2456200.5 epoch) provided by JPL's HORIZONS 
system (Giorgini et al. 1996; Standish 1998). In all the figures, t = 
coincides with the JD2456200.5 epoch. Relative errors in the total 
energy at the end of the simulations are < 5 x 10~^^. In addition 
to the calculations completed using the nominal orbital elements 
in Tables [T] and |2] we have performed 20 control simulations (for 
each object) using sets of orbital elements obtained from the nomi- 
nal ones within the accepted uncertainties (3a). The sources of the 
Heliocentric Keplerian osculating orbital elements and uncertain- 
ties are the JPL Small-Body Databas^EI and the AstDyS-2 portallfl 
The results of these control calculations are consistent with those 
from the nominal ones. In order to study the long term stability of 
these newly found Trojans we use the Regularized Mixed Variable 
Symplectic (RMVS) integrator (SWIFT-RMVS3) which is part of 
the SWIFT package (Levison & Duncan 1994). The SWIFT pack- 
age is available from H. Levison web siteQ The physical model 
here is similar to the previous one but the Earth and the Moon are 
replaced by the barycentre of the system. One group of calcula- 
tions does not include Mercury or the Pluto-Charon system, uses a 
timestep of 9 d and has been followed for ±4.5 Gyr with six con- 
trol orbits per object; the second group includes Mercury and the 
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Pluto-Charon system, uses a timestep of 3 d and runs for ±2.25 
Gyr. Relative errors in the total energy are now < 1.5 x 10~*. Re- 
sults from RMVS and Hermite are consistent within the applicable 
timeframe (±2 Myr). 



3 KNOWN MARS TROJANS 

In this section, we present the results of the short-term integrations 
of the orbits of the known Mars Trojans (nominal orbital elements 
in Table[T). These calculations are presented here because there are 
no references including a dynamical analysis of all of them; these 
results are also used to validate our simulations against previous 
work. Our results are consistent with those from other authors. 

5261 Eureka. Discovered by D. H. Levy and H. E. Holt on 
1990 June 20 during the course of the Mars and Earth-crossing 
Asteroid and Comet Survey conducted by E. M. and C. S. Shoe- 
maker (Holt et al. 1990), the first known images of this object date 
back to 1979 November and its dynamics has been studied in de- 
tail for more than 20 yr. Its diameter is 1.3-2.6 km (Trilling et al. 
2007; Burt 2012), its rotational period is 2.69 h and it has a sec- 
ondary companion with an orbital period of 16.93 h (Burt 2012). 
The period of its L5 Trojan motion varies between 1300 and 1400 
yr (Mikkola et al. 1994) as a result of the changes in inclination; the 
orbit is stabilized by the perturbation of the neighbouring planets. 
Its relative mean longitude is modulated with a major periodicity 
of nearly 200 000 yr, which is the circulation period for the longi- 
tude of perihelion, f2 ± a;. Our calculations confirm that it is an L5 
Mars Trojan as well as all the results in Mikkola et al. (1994), see 
Fig-Hi Its Ar oscillates around -60° with an amplitude of 11° and a 
libration period of 1365 yr. Its eccentricity oscillates mainly due to 
secular resonances with the Earth and the oscillation in inclination 
appears to be driven by secular resonances with Jupiter. 

(101429) 1998 VF31. This 800 m minor body (Trilling et al. 
2007) was nominated as L5 Trojan by Tabachnik & Evans (1999) 
shortly after its discovery. We also confirm that it is an L5 Mars 
Trojan. Its \r oscillates around -60° with an amplitude of 45° and 
a libration period of 1425 yr, see Fig. [T] The libration amplitude is 
much larger than that of Eureka and, consistently, its orbit appears 
to be less stable. As in the previous case, oscillations in eccentric- 
ity are driven by the Earth and those in inclination by Jupiter. Fig.[T] 
clearly shows that 101429 is not a dynamical analogue of Eureka; 
the evolution of its Ar, a and e is not similar to that of Eureka. 
The intrinsically different nature of these two objects has been fur- 
ther confirmed by spectroscopic results (Rivkin et al. 2007; Lim et 
al. 2011). Its origin points to a possible capture event 3.9 Gyr ago 
(Rivkin et al. 2007). 

(121514) 1999 UJ7. This object is first mentioned as an L4 
Mars Trojan by Rivkin et al. (2003); they also point out that its 
visible spectrum is different from those of Eureka or 101429. It is 
likely a P-asteroid and the largest of the known Mars Trojans with 
a rotational period of 21.2 h (Burt 2012). Our calculations confirm 
that this object is following a tadpole orbit east of Mars so it is an 
L4 Mars Trojan. Its A,, oscillates around -1-60° in very wide, asym- 
metric loops with an amplitude of 77°, the largest of all the objects 
discussed here, and it has a libration period of 1500 yr, see Fig.|2] 
Simulations (Mikkola et al. 1994) show that, for Mars, stable ex- 
cursions about L4 and L5 as large as 80° occur on time-scales of 
millions of years which is consistent with our own findings. Oscil- 
lations in both eccentricity and inclination are also observed. 

2001 DH47. We also confirm that this object is a robust L5 
Mars Trojan. Its Ar oscillates around -60° with an amplitude of 



Three new stable L5 Mars Trojans 3 



Table 1. Heliocentric Keplerian orbital elements of known Mars Trojans. Values include the 1-cr uncertainty. (Epoch = JD2456200.5, 2012-Sep-30.0; J2000.0 
ecliptic and equinox. Source: JPL Small-Body Database and AstDyS-2.) 







Eureka 


101429 


121514 


2001 DH47 


311999 


semimajor axis, a (au) 




1.523 536 2 


1.524 181 4 


1.524 424 5 


1.523 798 6 


1.523 756 9 






±0.000 000 005 


±0.000 000 02 


±0.000 000 12 


±0.000 000 2 


±0.000 000 03 


eccentricity, e 




0.064 767 3 


0.100 31 1 5 


0.039 184 3 


0.034 659 6 


0.053 978 1 






±0.000 000 08 


±0.000 000 13 


±0.000 000 12 


±0.000 000 6 


±0.000 000 6 


inclination, 2 (°) 




20.283 11 


31.298 20 


16.749 27 


24.398 79 


18.620 98 






±0.000 007 


±0.000 02 


±0.000 012 


±0.000 05 


±0.000 02 


longitude of ascending node, Q (° ) 




245.069 82 


221.331 79 


347.396 04 


147.432 50 


282.499 55 






±0.000 014 


±0.000 02 


±0.000 02 


±0.000 04 


±0.000 04 


argument of perihelion, uj (°) 




95.414 85 


310.573 87 


48.417 49 


17.615 


176.873 4 






±0.000 07 


±0.000 1 1 


±0.000 2 


±0.002 


±0.0002 


mean anomaly, IVI C^) 




236.388 24 


32.830 35 


321.501 44 


53.521 


120.741 26 






±0.000 07 


±0.000 13 


±0.000 2 


±0.003 


±0.000 14 


absolute magnitude, H 




16.1 


17.4 


16.9 


19.7 


17.8 




121514 (1999 UJ7) 




time (Myr) 

Figure 2. Same as Fig.[T]but for (121514) 1999 UJy, the only known L4 
Trojan. 



Table 2. Heliocentric Keplerian orbital elements of the new objects studied 
in this research. Values include the l-cr uncertainty when available. (Epoch 
= JD2456200.5, 2012-Sep-30.0; J2000.0 ecliptic and equinox. Source: JPL 
Small-Body Database and AstDyS-2.) 







2011 SCigi 


2011 SL25 


2011 UNg3 


a (au) 




1.523 809 1 


1.523 931 


1.523 782 






±0.000 000 5 




±0.000 002 


e 




0.044 121 2 


0.114 465 7 


0.064 627 






±0.000 001 3 




±0.000 007 






18.744 79 


21.497 41 


20.363 2 






±0.000 05 




±0.000 3 






5.799 81 


9.431 98 


223.573 






±0.000 03 




±0.000 2 


Ul (°) 




196.400 


53.252 46 


165.241 






±0.005 




±0.005 






4.422 


146.943 43 


192.529 






±0.006 




±0.004 


H 




19.3 


19.5 


19.9 



2001 DH47, see Fig.[T] Out of all Mars Trojans, it currently has the 
smallest relative (to Mars) semimajor axis, 5.9 x 10^^ au. 



11° and a libration period of 1365 yr. Together with Eureka, it has 
the smallest libration amplitude of all the objects studied here. Its 
orbital behaviour is also very similar to that of Eureka, see Fig.[T] 
(311999) 2007 NS2. Again, we face a robust L5 Mars Tro- 
jan. Its Xr oscillates around -60° with an amplitude of 14° and a 
libration period of 1 3 1 yr. Its orbit is similar to that of Eureka and 



4 NEW MARS TROJANS 

The three new Trojans follow tadpole orbits around Mars' L5 and, 
as in previous cases, they exhibit oscillations in both eccentricity 
and inclination, see Fig. [3] Their distances to the Earth are always 
> 0.4 au and their arguments of perihelion uj circulate. 
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Figure 3. Same as Figs.HandlSbut for 201 1 SC191, 201 1 SL25 and 201 1 UNea. 



2011 SC191. It was originally discovered on 2003 March 21 by 
the Near-Earth Asteroid Tracking (NEAT) project at Palomar and 
named 2003 GX20, then was lost but re-discovered on 201 1 Octo- 
ber 31 by the Mt. Lemmon Survey (Pettarin et al. 201 1). Its orbit is 
robust; 45 observations from the years 2003, 2009 and 201 1 with an 
arc length of 3146 d have been used to compute it. Its orbital incli- 
nation is significant (~ 19°) but its eccentricity is low (0.04); both 
values are similar to those of known Mars Trojans. It is a relatively 
bright object, suitable for spectroscopy, with a typical apparent vi- 
sual magnitude of nearly 20. Our calculations indicate that it is a 
stable L5 Mars Trojan, see Fig.|3] Its relative mean longitude oscil- 
lates around -60° with an amplitude of 18° and a libration period of 
1300 yr, see Fig.ll] These values are comparable to those of the first 
Mars Trojan, Eureka (see Section 3); the evolution of their orbits is 
also very similar, see Fig. [T] The orbit appears to be very stable on 
Myr time-scales. As in the case of Eureka, \r is modulated with a 
major periodicity of nearly 200000 yr. 

2011 SL25. It was discovered on 2011 September 21 at the 
Alianza S4 Observatory on Cerro Burek in Argentina. Its orbit has 
been computed using 76 observations with an arc length of 42 d; 
even if somewhat reliable, it is in need of further observations in 
order to make it as robust as that of 2011 SC191. Because of its 
larger eccentricity, it can get brighter than 201 1 SC191 at nearly 19 
mag. Our calculations indicate that it is a stable L5 Mars Trojan, 
see Fig. [3] Its \r oscillates around -60° with an amplitude of 18° 
and a libration period of 1400 yr, see Fig.|4] 

2011 UNea. It was originally discovered on 2009 September 
27 by the Mt. Lemmon Survey and named 2009 SA170. After being 
lost, it was re-discovered on 2011 October 21, again by the Mt. 
Lemmon Survey. Its orbit is well known; 64 observations with an 
arc length of 793 days. Its orbital behaviour is similar to that of 
2011 SC191. Our calculations also indicate that it is a stable L5 
Mars Trojan, perhaps the most stable of this set of three, see Fig. 
[3] Its Ar oscillates around -60° with an amplitude of 14° and a 
libration period of 1350 yr, see Fig. |4] which are comparable to 
those of Eureka (see Section 3). 



5 DISCUSSION AND CONCLUSIONS 

The simulated time span studied so far (4 Myr) is not long enough 
to answer the astrophysically critical question of whether the newly 
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Figure 4. Left-hand panels: details of the evolution of the relative mean 
longitude for the three objects presented here in the time interval (-10 000, 
10000) yr. Right-hand panels: one tadpole loop starting at t = yr in the 
coordinate system corotating with Mars; tadpole loops are the superposition 
of multiple epicyclic loops. 



found Trojans may have been trapped at the Lagrangian L5 point 
of Mars since the formation of the Solar system. To investigate the 
long-term stability of the three objects we use the SWIFT-RMVS3 
integrator. Our ±4.5 Gyr integrations indicate that the lifetime of 
the Trojan orbits of the new objects is equivalent to or longer than 
the Solar system age (see Fig.|5}. They all could be primordial ob- 
jects. We also confirm the long-term stability of the four Mars Tro- 
jans studied by Scholl et al. (2005) and conclude that 311999 is 
also stable on Gyr time-scales. However, our calculations (includ- 
ing control orbits) suggest that 101429 and 121514 are not primor- 
dial but were captured about 4 Gyr ago. Some of the control orbits 
associated with 2011 SL25 are not stable over the entire simulated 
time. Our ±2.25 Gyr integrations show stability for all the eight 
objects. 

In this Letter, we have presented solid dynamical evidence in 
the form of A^-body numerical integrations for three new stable 
Mars Trojans and also provided a comprehensive summary of the 
dynamics of the previously known objects. After this study, the list 
of robust Mars Trojans includes eight objects. Unfortunately and 
with the exception of the first three objects, there is no consen- 
sus among Solar system researchers on the robust Trojan nature of 
2001 DH47 and 31 1999 in spite of both having reliable orbits. Even 
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Figure 5. The long-term evolution of the libration amplitude of all known 
Mars Trojans. All the objects appear to be stable on Gyr time-scales. 



the MFC still lists just three objects as genuine Mars Trojans0The 
A'^-body calculations completed for this research clearly show that 
the eight objects are long-term stable, at least for several Gyr. We 
also conclude that 101429 is by no means a dynamical analogue of 
Eureka - their orbital evolution being too different - which is con- 
sistent with the rather different spectroscopic signature found by 
Rivkin et al. (2007). Two of the three new L5 Trojans move in or- 
bits very similar to that of Eureka which, in the absence of surface 
composition data, may point to them (and several of the previously 
known Trojans) being the by-products of the disruption of a larger 
body. Eureka and 101429 could be coUisional fragments of larger 
bodies as they are highly differentiated (Rivkin et al. 2007). The 
binarity and relatively fast rotation of Eureka (Burt 2012) as well 
as the apparently strong asymmetry between the size of the L4 and 
L5 populations (one versus seven) also point in that direction. Al- 
though our calculations indicate that most of the currently known 
and new Mars Trojans are dynamical analogues of Eureka, spec- 
troscopy should settle the question of them also having a common 
chemical signature with Eureka, an R chondrite (Lim et al. 2011). 
If, as suggested, they are fragments of a large body, further ob- 
jects with similar properties should be found (Marzari et al. 1997). 
Strategies to search for additional Mars Trojans have been recently 
discussed by Todd et al. (2012). 
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NOTE ADDED IN PRESS 

After this work was accepted by MNRAS Letters, a relevant paper 
was submitted by A. A. Christou to astro-ph (astro-ph/1303.0420). 
He arrives to similar conclusions regarding the three objects dis- 
cussed here but using different techniques (Christou 2013). 
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